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ABSTRACT 

Aims. Chromatic effects are usually associated with refractive optics, so reflective telescopes are assumed to be free from them. We show 
that all-reflective optics still bears significant levels of such perturbations, which is especially critical to modern micro-arcsecond astrometric 
experiments. 

Methods. We analyze the image formation and measurement process to derive a precise definition of the chromatic variation of the image 
position, and we evaluate the key aspects of optical design with respect to chromaticity. 

Results. The fundamental requirement related to chromaticity is the symmetry of the optical design and of the wavefront errors. Finally, we 
address some optical engineering issues, such as manufacturing and alignment, providing recommendations to minimize the degradation that 
chromaticity introduces into astrometry. 

Key words. Astrometry, Methods: data analysis, Space vehicles: instruments, Techniques: high angular resolution, Telescopes. 



1. Introduction 

Chromatism is usually defined for refractive optics as an aber- 
ration due to the light dispersion of the glass with the refractive 
index; this induces a perturbation on the image profile, depend- 
ing upon the source color. In astrometry, it provides a variation 
in the apparent star position, i.e. an astrometric error. 
This effect can be reduced by using more complex refractive 
systems (doublets, triplets, etc.), taking advantage of the differ- 
ent dispersion of glasses to achieve a certain degree of compen- 
sation. For most purposes, the only way to avoid chromatism is 
considered to be adoption of a completely reflecting design. 
Here, we show that, at the demanding level of modern astrom- 
etry, reflective optics is still affected by significant chromatic 
effects, which we will refer to as "chromaticity. We investi- 
gate the source of the chromaticity and provide recommenda- 
tions for its minimization in optical design, manufacturing, and 
alignment. 

Currently, space astrometry experiments are being designed 
and implemented with the goal of micro-arcsecond (hereafter, 
yuas) measurements. We will refer to the framework of the Gaia 
mission, approved within the space science program of ESA 
for launch before 2012. Built upon the implementation of the 
Hipparcos (ESA 1997 1 concept, with the benefit of modern 
technology and more advanced astrophysical understanding, 
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Gaia aims to measure absolute position, parallax and annual 
proper motion of ~ 10 9 objects with a typical accuracy of ~ 
15 //as for V=15 stars, with survey completeness to V=18 and 
limiting magnitude V=20 (Perrymanet al. 2001). Hipparcos's 
payload was already affected by chromaticity at the mas-level, 
and during Gaia's design phase it soon became apparent that 
its all-reflective optics would also be affected by chromaticity 
at a similar level ( Latt anzi et al. 19 98). Therefore chromaticity 
in Gaia must be suppressed by more than two orders of magni- 
tude. 

We describe how the aberrations of a realistic all-reflective op- 
tical system can lead to chromaticity as a consequence of light 
diffraction, even when propagation occurs in a practically non- 
dispersive medium as would be the case for payloads sent out- 
side the Earth's atmosphere. Although our discussion focuses 
on Gaia, the concepts and considerations developed here about 
sources of chromatic errors and how to deal with them can ap- 
ply to all experiments wishing to reach astrometric accuracy 
down to the mas-level and beyond. 

The aspects evidenced by our study of chromaticity are of crit- 
ical importance for compiling detailed error budgets and de- 
ciding on actions for controlling systematic errors. Also, the 
degree of detail in instrument and measurement models re- 
quired to keep chromaticity at the //as-level can help prevent 
(during design), monitor (throughout operations), and correct 
(in data reduction) other possible sources of systematic errors 
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Gaia-like optical layout (Baseline Configuration) Positions: 1-2 DS 28-Apr-05 Gaia-like optical layout (Baseline Configuration) Scale: 0.09 DS 24-Feb-05 

Fig. 1. A Gaia-like astrometric payload layout. Left: the set of two telescopes combined at mirror M4 (planar drawing) Right: 
individual telescope layout. LOS A and LOS B are the two lines of sight implemented in this two-telescope configuration. 



associated with an image profile and its temporal and spatial 
variations. A possible method for chromaticity correction is de- 
scribed in Gai & Cancelliere (2005). 

In Sect. 2, we recall the relevant basic principles of a global 
astrometric mission like Gaia 1 . In Sect. 3, we introduce chro- 
maticity for an all-reflecting optical system through a simpli- 
fied analytical model, which is also utilized to prove the linear 
dependence of chromaticity on wavelength. Section 4 presents 
the chromaticity analysis for a realistic representation of the 
Gaia optical system. Section 5 provides prescriptions for min- 
imizing chromaticity in optical designs. In Sect. 6, we present 
the field-of-view (FOV) distribution of chromaticity and me- 
thods for its reduction. Finally, conclusions are drawn in Sect.7. 



2. Telescopes for global astrometry 

Global astrometry is the direct measurement of angular po- 
sition, parallax, and proper motion, i.e. of those fundamen- 
tal parameters that define, along with the line-of-sight (LOS) 
velocity, the location and status of motion of every object 
in space. The success of the Hipparcos mission established 
that the only meaningful way to achieve global astrometry 
is from a space-born platform. Suitable observing strategies 
and a sufficiently long operational life-time would allow the 
whole sky to be covered repeatedly thus providing the means 
to derive positions and their variations with time. As the re- 
sult of renewed interest in global astrometry by the astronom- 
ical community worldwide, several missions are being con- 
sidered by the major space astronomy communities. In Japan 
studies are underway on the near-infrared astrometry mission 
JASMINE (IGouda et al. 2005> for a possible launch in 2014. In 
the United States, the NASA funded concept study for the pos- 
sible development of the Origins Billion Star Survey (OBSS) 
is close to completion (John ston et al. 2 005 1, while teams at Jet 
Propulsion Laboratory (JPL) are well into the development of 
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SIM (Space Interferometry Mission), due for launch in 2011 
dMarr2003> . 

In Europe, following formal approval in 2000, the European 
Space Agency (ESA) is developing its new astrometric mission 
Gaia, which is based on the same two-FOV's scanning satellite 
concept implemented on Hipparcos, and is meant to increase 
the measurement capabilities of its predecessor by about two 
orders of magnitude. 

Two different optical systems were considered during the 
study phase, so that we will not know the optical configura- 
tion that will fly on Gaia before final selection of the indu- 
strial contractor (expected by early 2006). However, this is not 
a problem for the scope of this work as both designs are based 
on all-reflective optics, and we can still refer to one of them 
without loss of generality. We chose to refer to the configu- 
ration known in Gaia jargon as the "Baseline Configuration" 
(Perry man et al. 2001). It consists of the combination of two 
off-axis, three-mirror monolithic telescopes mounted on a com- 
mon optical bench, as shown in Fig. 1 (left). The individual 
telescope has an effective focal length (EFL) of 46.7m and a 
rectangular collecting area of 1.4x0.5 m 2 , resulting in a central 
lobe size (analogous to the Airy diameter for circular pupils) of 
~ 250 x 700 mas at the effective wavelength A e g = 700 nm. We 
constructed a computer model of the Baseline Configuration by 
utilizing the CODE V optical modelling package (ORA 1995 1, 
and optimized the system over a FOV of 0°.66 x 0°.66 by 
closely following the procedure described in Loreggia et al. 
(2004). With reference to Fig.^right), the optical train of each 
telescope consists of seven elements: three mirrors (Ml, M2, 
M3) with optical power, three flat mirrors (M4, M5, M6) for 
beam-folding, and the focal plane (FP). The FOV covers the 
angular ranges [-0°.33, 0°.33] in the along-scan direction and 
[0°.20, 0°.86] in the across-scan direction; the off-axis design is 
required to avoid vignetting. One of the possible implementa- 
tions investigated for the Gaia FP assembly calls for an array of 
17 (along-scan) x 10 (across-scan) CCDs working in time de- 
lay integration (TDI) mode ( Sai nt-Pe et al. 2000 1. Each along- 
scan strip of devices hosts: two chips for the functions of tar- 
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get detection and confirmation, eleven CCDs for the repeated 
astrometric measurement, four units for broadband photome- 
try. The rationale for photometry is both astrophysical (stellar 
classification) and technical (chromaticity correction discussed 
below). 

During each transit, the image of a target on any of the eleven 
astrometric CCDs is binned across scan, and centered along 
scan, providing independent one-dimensional location esti- 
mates of comparable precision (IGai et al. 20011 and referen- 
ces therein). The composition of repeated along-scan measure- 
ments, taken for every location on the sky with different orien- 
tations, allows the determination of the bi-dimensional angular 
coordinates in a common reference system ("global sphere re- 
construction"). The evolution of the apparent source position 
during the mission lifetime (~ 5 years) yields absolute parallax 
and proper motion. 

Note that the relative orientation of the two telescopes, i.e. the 
base angle (BA), is a key factor of the measurement process: 
the on-sky separation between two stars entering the two dif- 
ferent FOV's is obtained from the measured difference in FP 
coordinates plus the BA; therefore, configuration perturbations 
of a few nm, negligible as to the optical response, are critical 
to astrometry as they can induce errors up to several tens //as. 
Laser-based metrology lines are expected to monitor and (if 
necessary) to control the BA stability that also sets stringent 
constraints on pay load thermal stability (Ga rdiol et al. 2004 1. 

3. Chromaticity 

3.1. A basic analytical model 

Chromaticity in all-reflective optical systems can be defined as 
the shift in the FP location of the images of two sources of dif- 
ferent spectral types sharing the same direction in object space. 
Also, chromaticity is an intrinsic characteristic of diffraction 
that cannot be neglected when the targeted accuracy level is at 
the mas level and higher, as was the case for Hipparcos and will 
be for Gaia. 

We verified the presence of chromaticity is mainly related 
to wavelength; sampling of the telescope point spread function 
(PSF) and dimensions of the read-out region of the FP play a 
secondary role in the actual amount of chromaticity displayed 
by the realization of any all-reflective optical design. 
Also, we will see that for a symmetric PSF there is no chro- 
maticity apart from possible residual effects due to sampling 
and other noise sources (e.g. photon noise). On the other hand, 
for an asymmetric 2 PSF, most of the chromatic effect, A, can be 
modelled for every location estimators tested as a linear func- 
tion of the change in effective wavelength 5/t e ff, characteriz- 



2 In quite different optical configurations, the location of the peak of 
the diffraction pattern appears to be stable with respect to the spectral 
distribution of the celestial objects; this is consistent with the fact that, 
in case of reasonably limited aberrations, this maximum is coincident 
with the chief ray of geometric optics, which by definition is insen- 
sitive to wavelength variations. Therefore, the diffraction peak would 
appear as a good candidate for the PSF "photo-center" (location), if 
only the PSF were sampled with infinite resolution. In practice, it is 
impossible to measure it directly because of the finite pixel size. 



ing the spectral difference of the sources, through a coefficient 
g(x,y) that depends on the local shape of the PSF, f{x,y), i.e. 
on the contributions of the aberrations over the FOV: 



A = g(x,y) ■ 6Aes, 



(1) 



where (x,y) are linear coordinates on the FP. 

We illustrate relation Q with the help of ad hoc analy- 
tical model. We choose an asymmetric, one-dimensional bell- 
shaped curve, in the form of a deformed Gaussian profile, i.e. 
of different width on the opposite sides of its peak value (at 
0) 3 : 
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e being the "asymmetry" factor (e << <xq). For continuity of 
the PSF, we must set 



cr(0) = o~q 



/(0) 



Ao~o V27T 



(3) 



Note that the position of the maximum is the origin, indepen- 
dent of A. 

For two different wavelengths A\ < A-i, the normalization con- 
dition requires that the maximum value decreases as the curve 
width increases. This agrees with the physical condition that 
in the diffraction limit, the PSF of an optical system increases 
in linear size with the wavelength. In particular, the PSF is a 
function of the wavelength A, the position on the FP y, and of 
the characteristic linear dimension D of the optical system, in 
a well-defined form required for dimensional reasons: 



PSF(y,A,D) = PSF 



(?)■ 



(4) 



This is the typical form of the argument of the diffraction pat- 
tern for the analytical reference cases solved in the literature. 
Therefore, a change in wavelength corresponds, for the PSF, to 
a scaling of the FP coordinates; any parameter weighted with 
this distribution will be affected by such variation. 
Let us evaluate the behavior of, e.g. the center of gravity (COG) 
estimator, for the model in Eq. over the region of interest 
[—a, a]: 



yf(y) dy = £ 

a 



O~0 



4e + ^Jc£ll . exp • 



.1 r s i 2 

2 \_Mcr -e)\ 



(5) 



where the approximation is valid for a > ctqA. The COG rep- 
resents one of the possible estimators of the photo-center loca- 
tion. In the Gaia literature, for example, the photo-center loca- 
tion of a realistic, digitized image is named "centroid" for each 
estimator used. Nevertheless, in this paper we retain the expres- 
sion COG for ease of the reader. The photo-center is linearly 
dependent on the selected wavelength through an asymmetry 



3 As we make reference to a Gaia-like system in this work, we con- 
sider the y axis for consistency with most of the literature on Gaia. 
There the y-axis is the axis of the one-dimensional location process 
and the axis along which the scanning occurs. 
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factor e, representing the difference between the two sides of 
the PSF, and the chromaticity relation Q becomes: 



A = 6<y) = 



4e 
V2^ 



SX 



eff- 



(6) 



In a more general case, the parameter e can be considered 
a slowly varying function of the focal plane coordinates, i.e. 
e(x,y), to account for the variation of the optical response 
(aberration distribution) over the FOV. Also, different location 
estimators would generate functional forms different from that 
in Eq. l|6}. 

3.2. Asymmetric PSF with coma 

The results of the previous section can be verified, in a more 
general way, from the Fraunhofer integral for a circular aper- 
ture of radius a: 



U(u, v) = C 



o Jo 



exp[i(k<&-vpcos(6->p)+up 2 /2)]pdpd6, (7) 



where C is a normalization constant, (p, 9, z) and (r, i/r, z) are the 
cylindrical coordinates in the exit pupil and focal plane (image) 
spaces, respectively. Then O is the term of aberration, < p < 
1, and the dimensionless variables u and v are defined as: 



2n a 



2n a 

v — r 

A R 



— -yjx 2 +y 2 , 



(8) 



R being the radius of the Gaussian sphere (=EFL). Following 
the treatment in Born and Wolf (1980), we assume working 
with no piston (u = 0) and with the tangential coma as the 
only aberration term; this type of aberration is expressed in the 
Nijboer-Zernike representation as $ = A c (3p 3 - 2p) cos(6>) (A c 
is a coefficient that quantifies the amount of coma in wave- 
lengths). In this case, it can be shown, that the solution of Eq. 
Q takes the form of a series: 



U(u, v, i/O = t/(0, v, ip) = C[U (0, v, ifr)+ 

+ (iaro3i)t/i(0,v,t/0 + (/a-o3i) 2 f/2(0,v,i/0 + ...], 



(9) 



where ao3i = 2ttA c /A and the functions t/j are given as a com- 
binations of Bessel functions of different orders, J^: 



t/o(0,v,«A) = 



2Ji(v) 



t/i(0,v,i/0 = icos(^) 



2J 4 (v) 



U 2 (Q,v,M = — 
2v 



\ji(y) - ^3(v) + ij 5 (v) 



=- T^Mv) + cos(2.A) 
2v 20 



(10) 



We used Eqs.© and J10I and the value 4 A c = 0.3 (in 
wavelengths) for the coma coefficient to compute the PSF. 
Hereafter we use the one-dimensional version, the Line Spread 

4 In practice, the optical system would be optimized to reduce 
coma. Consequently, the actual variation of the COG would, of course, 
be smaller. 




Fig. 2. The LSF obtained from the Fraunhofer integral with the 
coma term <L> (see text). 

Function (LSF), obtained from the PSF by across-scan integra- 
tion. The resulting curve is shown in Fig|5]for the two refer- 
ence spectral types considered, B3V and M8V in the case of 
a Gaia-like configuration with a = 0.7 m and R = EFL = 
46.7 m. Note that the intrinsic effective wavelengths of B3V 
and M8V stars are about 150 nm and 1300 nm, respectively. 
After taking realistic telescope transmission and detector re- 
sponse ( |Short & de Bruijne~2 003 ) into account, the effective 
wavelengths become Abw = 628 nm and Amhv = 756 nm. 
Using the discrete representation of the COG: 



JCOG - T: y^^nfn, 



(ID 



with f„ the signal intensity recorded by the pixel in y n , and 
F = Ytfn the total intensity, we verified that the variation of 
the COG position is linear with wavelength in the range of 
interest following the relation dycoc = -1.0195/1 + 0.6415, 
in agreement to the result derived in Sect. 13.11 The calcula- 
tion was carried on by making reference to a situation using 
LSF sampling and focal plane detectors, very similar to one 
of the options considered for Gaia tShort 201751 : single CCD 
units of about 50 mm in the scan direction, and pixels of 10 pm 
in the same direction. The difference between the COG loca- 
tions of the LSFs of our two reference star types amounts to 
|COG B3 v - COGmsvI = 0.121mas. 

With our PSF modelling and measuring code we were able to 
investigate the effects of limited detector dimensions (read-out 
region coincident with the detector dimension) and sampling 
resolution (pixel size). We went on to extend the detector di- 
mension and to shrink the sampling step by two orders of mag- 
nitude (CCD units of 50 cm on a side with 0.1 pm pixels) 
and repeated the COG location measurements, which yielded: 
|COGb3v - COGmsvI = 0.122 mas. This is a variation of only 
1% compared to the value above. 

These experiments suggest, within the limits probed (the di- 
mension of the "unlikely" CCD simulated is very large but not 
infinite, and its pixels small but not infinitesimal), that chro- 
maticity is an intrinsic property of all-reflective optical systems 
and that it can be approximated by a linear function of source 
effective wavelength. Therefore whatever the detector geom- 
etry and its spatial extension, different spectral type stars, set 
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in the same nominal position in object space, do not have the 
same estimated position on the FP, even with the simplest pos- 
sible algorithm; i.e. the same location on the sky is not uniquely 
mapped on the focal plane. 

Chromaticity, a color-dependent position variation, is an impor- 
tant contribution to the systematic error and must be removed 
to preserve the desired mission accuracy. 

4. Chromatic astrometric error 

Using our CODE V representation of the Gaia Baseline 
Configuration, we could derive Gaia-like images of the two 
B3V and M8V sources (all aberrations terms are now taken 
into account), for any representative off-axis field position. The 
sources are always intended to be in the same location in ob- 
ject space, i.e. along the same direction on the sky. An example 
of the one-dimensional optical signals are shown in the upper 
panel of Fig.|5J normalized to the peak value; it is apparent that 
the M8V image (dashed line) is not just a geometrically scaled 
version of the B3V case (solid line). The different spectral con- 
tent weighs the design aberrations (i.e. the residual aberrations 
after final optimization of the optical design), and the diffrac- 
tion image is also affected by profile variations, in addition to 
the magnification associated to the intrinsic diffraction factor 
A/D. The signal difference is shown in Fig. [5] (lower): it is 
not symmetric, as would be the case for perfect scaling, and 
reaches 15% of the initial peak value. For consistency with the 
treatment of the previous section, the coordinate system is de- 
fined so that the COG location of an ideal, "non-aberrated", 
PSF has abscissa equal to 0; the large common mode displace- 
ment (the PSF's do not peak at 0) is mostly due to classical 
distortion. It can be shown that different algorithms are af- 
fected by similar, although not coincident, chromatic effects. 
Nevertheless, because of the effective image profile of the Gaia 
telescope, the difference in performance is relatively small (few 
percent), so that we continue to prefer the simplicity of the 
COG algorithm throughout this work. 

The application of the COG method for centering the Gaia-like 
PSF's of the two reference stars yielded 20. 112 fim for the blue 
star and 20.055 pm for the red star, respectively, a difference of 
57 nm. At the optical scale of Gaia, about 4"/mm, this results 
in more than 200 yt/as, i.e. comparable with the random location 
error for stars brighter than V = 15 mag. 

4.1. Monochromatic and polychromatic PSF's 

It is possible to simplify the formalism used in sect. I3.2l bv re- 
calling that the PSF can be derived by taking the square modu- 
lus of the complex amplitude response function (ARF), which 
in turn is derived as the Fourier transform of the generalized 
pupil function P s (Goodman 1996), i.e. 



ARF(x,y; A) = q(A) X F{P g (p, 9)} , 



(12) 



where (x, y) are coordinates on the focal plane and (p, 9) are 
taken on the pupil. 

The Fourier integral in Eq. fl!2i is performed over the rect- 
angular pupil of the Gaia telescope, where the factor q(A) is 




-20 -10 10 20 30 40 
Along scan position \jum] 



50 60 



Fig. 3. Upper: PSF for B3 V (solid line) and M8 V (dashed line) 
stars in a representative off-axis field position; lower: PSF dif- 
ference. 

constant. The pupil function P g is expressed in terms of the 
wavefront error (WFE) function, W(p, 9), via the relation 



P g (p, 9; A) = exp 



In 

i— W(p,0) 
A 



(13) 



Note that Eq. dl 31 clearly shows how at longer wavelengths 
the effect of the WFE is reduced, as the phase contribution is a 
smaller fraction of the period. 

The expression for the monochromatic PSF, PSF m , is therefore 



PSF m (x,y,A) = \ARF(x,y;A)\ 2 



(14) 



The local polychromatic PS F(x, y) for a star with a given spec- 
tral distribution is 



PSF(x,y) 



!' 



u)(A)xPSF m (x,y,A)dA 



(15) 



The weight function a>(A) combines the source spectral distri- 
bution, the instrument transmission, and the quantum efficiency 
of the detector; then the polychromatic PSF can be normalized 
to the expected average number of detected photons. 
CODE V can compute the WFE distribution associated to a 
given field position for a selected optical configuration; it is 
thus possible to build the PSF for any desired source, using Eqs. 
H2\ - H5\ . We use a numeric implementation of this model 
( Bus onero et al. 200 5 1 to provide the description of the imag- 
ing performance of the Gaia telescope with all the aberration 
terms given by CODE V. Additional contributions describing 
the realistic detector response and the effect of TDI observa- 
tions can be included, as well as the across-scan binning used 
for Gaia, to build the recorded signal. For any source of known 
spectrum, it is possible to derive its effective wavelength, the 
detected signal, and the chromaticity with respect to the se- 
lected reference spectral type using the above expressions. In 
this way, the measurement of Gaia can, in principle, be made 
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Table 1. Effects of individual Standard Zernike terms vs. the non-aberrated (ideal) case. 



Term 


Standard 


RMS WFE 


Image RMS 


COG displacement 


Chromaticity 




Zernike 


[nm] 


width increase [%] 


[mas] 


[pas] 


no. aber. 


/ 














1 


1 














2 


pcos(S) 


11.37 











3 


p sin(#) 


32.44 


0.1346 


34.29 


97.95 


4 


p 2 cos(20) 


13.02 


1.263 


0.030 


-2.22 


5 


2p 2 - 1 


26.03 


4.927 


0.119 


-9.27 


6 


p 2 sin(20) 


26.48 


0.3199 


21.35 


-359.7 


7 


p 3 cos(36>) 


11.76 


1.881 


0.045 


-3.651 


8 


(3p 3 - 2p) cos(0) 


12.51 


1.881 


0.045 


-3.691 


9 


(3p 3 - 2p) sin(0) 


32.49 


1.455 


-33.38 


779.6 


10 


p 3 sin(3(9) 


13.61 


0.829 


14.34 


-118.3 


11 


p 4 cos(40) 


8.928 


1.079 


0.028 


-1.042 


12 


(4p 4 - 3p 2 )cos(20) 


23.15 


5.026 


0.144 


11.64 


13 


(6p 4 _ 6p 2 + 1) 


34.21 


9.965 


0.262 


20.22 


14 


(4p 4 -3p 2 )sin(26») 


38.69 


2.083 


-36.64 


1019 


15 


p 4 sin(4#) 


6.423 


-0.008 


0.125 


-652 


16 


p 5 cos(50) 


5.597 


0.3312 


0.097 


0.6501 


17 


(5p 5 -4p 3 )cos(36») 


246.53 


7.817 


0.234 


26.46 


18 


(10p 5 - 12p 3 + 3p)cos((9) 


25.99 


6.186 


0.175 


14.92 


19 


(10p 5 - 12p 3 + 3p)sin(6>) 


20.61 


3.329 


2.229 


-2142 


20 


(5p 5 -4p 3 )sin(36>) 


26.35 


4.571 


-18.66 


1372 


21 


p 5 sin(50) 


5.376 


-0.169 


-1.346 


-393.4 



achromatic, as the position estimate is no longer function of the 
source spectral distribution. This correction requires both color 
information and knowledge of the local instrument response. 

4.2. Identification of the critical WFE terms 

The source of chromatic errors is the WFE, which is present in 
any nominal optical design with a finite FOV; manufacturing 
and alignment may only aggravate the problem. Besides, not 
every component of WFE contributes to chromaticity: below, 
we analyze the effect of individual aberrations and the related 
symmetry properties. 

We then investigate how to minimize, by design, the chromatic 
errors in any given field position, analyzing the effect of partial 
or total suppression of selected aberration terms. Finally, we 
proceed to evaluate the overall field properties of chromatic- 
ity, verifying the possibility of compensating for it by taking 
measurements over a complete crossing of the focal plane, a 
possibility offered by scanning instruments like Gaia. 
In optical engineering, typical expansions of the WFE are 
in terms of Zernike and Fringe Zernike polynomials, or- 
thogonal, and normalized functions, if mapped on a cir- 
cular pupil (Bor n & Wolf 19801 . In our analysis we have 
used expansions with 21 Zernike terms; the set of standard 
Zernike functions are listed in Table ^ (f° r Fringe Zernike 
equivalents see Born & Wolf 1980 1. Note that the representa- 
tion in Zernike polynomials is not optimal for non-circular 
pupils and more convenient expansions have been investigated 
JGardiol et al. 20051 

In the following, we show how to identify the main individ- 
ual contributors to chromaticity in the case of a realistic rep- 
resentation of the Gaia telescope regardless of the choice of 



the Zernike set, and to search for possible correlations between 
chromaticity and parameters like RMS WFE and image RMS 
width. 5 

Chromaticity is again evaluated by taking the difference of 
the COG positions obtained for the adopted reference spec- 
tral types B3V and M8V. We used both monochromatic PSFs 
at the reference wavelengths associated with the two spectral 
types and the full polychromatic representation of Eq. d!5l >. as- 
suming the two stars behave like blackbodies. The latter de- 
scription is more representative from the astrophysical stand- 
point, whereas the former is computationally much simpler. 
Preliminary experiments proved that the two methods can be 
considered equivalent when, as is the case here, the primary 
goal of the investigation is to evaluate of the properties of chro- 
maticity and not its precise calculation. 

The image RMS width is derived from the monochromatic 
PSF at A — 700 nm, roughly representative of solar type ob- 
jects and sort of midway between B3V and M8V types. The 
image COG displacement is always referred to the ideal, non- 
aberrated case, which has an image RMS width of 8.6425 fim. 
Each Zernike term is individually evaluated, with the coeffi- 
cient set to 0.1 (i.e. the small aberrations regime is assumed): 
then the WFE is built from the selected term, and the PSF com- 
puted accordingly to the above model. The results for the first 
21 terms are listed in Tabled 

In Fig.0] we show the impact of each individual aberration on 



5 The root mean square wavefront error, RMS WFE, is simply the 
square root of the quantity j WFE 2 (p,6) taken over the pupil. The 
image RMS width for, e.g. the y axis on the focal plane is computed as 

J f(y— < y >) 2 L(y)dy, where L{y) = J PS F(x, y)dx is the line spread 
function. 
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Fig. 4. Image COG variation (upper panel) and chromaticity 
variation (lower panel) vs. Standard Zernike aberration terms. 

the COG (i.e. an effect corresponding to classical distortion) 
and on chromaticity. Some terms strongly contribute to both, 
but in general there is little correlation. Similarly, we found that 
there is no simple relationship between RMS WFE, or image 
RMS width (both including the contribution of all aberrations) 
and chromaticity. 

The complex relationship between WFE, image RMS width, 
and chromaticity can be understood in terms of symmetry, 
since the aberrations are in general bidimensional, whereas the 
parameters relevant to the Gaia measurement are mostly one- 
dimensional, i.e. referred to the scan (y-axis) direction. Thus, 
specific aberration terms may contribute to the WFE signif- 
icantly, to across-scan image width and across-scan (x-axis) 
centroid displacement, with little impact on the location, noise 
and chromaticity, due to the across-scan binning of the images 
(see Sect.|2j. 

There is a trend toward increasing image RMS width with RMS 
WFE, but this is not a strict relation, and several cases of large 
WFE and small image width degradation are seen. This could 
be associated to aberrations inducing significant image degra- 
dation in the across scan direction only. The COG, in some 
cases, is affected by a large displacement associated to low 
chromaticity values. Thus, the images are translated with re- 
spect to the non-aberrated position by an amount that does not 
depend on the spectral distribution of the source. The case of 
classical distortion fits this description. 

From Fig. |4] we see that significant contributions to chro- 
maticity, i.e. > 100 //as, come only from aberrations (no. 3, 6, 
etc.) with a specific functional form: they are all the odd (si- 
nusoidal) functions of the angular coordinate, so that they have 
odd symmetry on the pupil plane. All of the even (cosinusoidal) 
terms do not provide net chromaticity; values below ~ 30 yuas 
are the result of limited precision in the calculation scheme that 
is implemented. 

A similar analysis was carried on with the Fringe Zernike poly- 



nomials. Again, the critical terms for chromaticity we identified 
to be those associated to the odd parity portion of the WFE, i.e. 
the sinusoidal terms, as for the Standard Zernike expansion. 
Therefore, in the case of Gaia, the terms that should be mini- 
mized, by design, manufacturing, and alignment are those with 
odd parity (anti-symmetric) relative to the across-scan (x) axis. 

5. Compensation by aberration selection 

We now turn our attention to verifying chromaticity with many 
aberration terms present. The analysis refers to FP position 
F4 (see Table EJl with coordinates 0°.20 (along-scan) x0°.33 
(across-scan), which is affected by a significant value of chro- 
maticity 6 . Six cases are considered: 

1 . non-aberrated (ideal) PSF, 

2. nominal case (from optical design), 

3. removal of all symmetric terms, 

4. removal of all anti-symmetric terms, 

5. removal of anti-symmetric terms AND scaling of symmet- 
ric terms to equivalent RMS WFE, 

6. subset of anti-symmetric terms 9 and 10. 

The corresponding values of RMS WFE, image RMS width, 
and chromaticity are listed in Table [3] The resulting PSF's 
(monochromatic at the reference wavelength of 700 nm) for 
four of the cases above are shown in Fig. [5] The nominal case 
provides an image RMS width that is reasonably close to the 
diffraction limit in spite of non negligible RMS WFE (40 nm, 
i.e. A/15 at A = 600 nm); the chromaticity is about 1 mas. 

Suppression of the symmetric aberrations provides some 
improvement to the RMS WFE. The variation of the image 
width is marginal and the initial chromaticity is mostly re- 
tained. Conversely, removing the anti-symmetric aberrations 
completely (case 4), we achieve significant improvement, as 
in the previous case, on the RMS WFE, but not nearly as much 
on the image RMS width; chromaticity, however, is reduced to 
zero, according to expectations. The RMS sum of WFE in cases 
3 and 4 restores the nominal value: the mutual orthogonality of 
symmetric and anti-symmetric function sets is preserved, even 
if this is no longer true for the individual functions within each 
set. Even when the symmetric aberrations are scaled to restore 
the initial WFE level of 40 nm (case 5), the chromaticity is still 
zero. 

Besides, when a random subset of anti-symmetric aberrations 
(9 and 10) is retained, together with the symmetric terms, with 
the nominal coefficients (case 6), the result has a value that is 
comparable to the chromaticity in the nominal case (case 2) but 
with the opposite sign. Futhermore it leads to a dramatic degra- 
dation of both WFE and image quality. This shows that the op- 
timization procedure in the ray tracing code actually achieves 
some partial compensation among different aberrations, pro- 
viding some, although limited, benefit to chromaticity. On the 
other hand, standard optical design optimization procedures 



6 As for the previous section, chromaticity is evaluated as the differ- 
ence between the monochromatic PSF at 628 nm and 756 nm, respec- 
tively. The chromaticity derived with the fully polychromatic model is 
about a factor 2 larger. 
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Table 2. RMS WFE, image RMS width (for the two reference spectral types), and chromaticity at different field angle locations 
over the focal plane. 





X position 
[degrees] 


Y position 
[degrees] 


RMS WFE 
[nm] 


Image RMS width 
for a B3V star [mas] 


Image RMS width 
for a M8V star [mas] 


Chromaticity 
[//as] 


Fl 


0.53 


0.00 


22.03 


39.27 


42.85 


0.00 


F2 


0.86 


0.33 


208.73 


41.80 


43.89 


-3071.23 


F3 


0.86 


-0.33 


208.73 


41.81 


43.89 


3070.16 


F4 


0.20 


0.33 


39.77 


41.93 


44.92 


2077.02 


F5 


0.20 


-0.33 


39.77 


41.93 


44.92 


-2077.02 


F6 


0.70 


0.17 


24.32 


38.51 


42.54 


124.94 


F7 


0.70 


-0.17 


24.32 


38.51 


42.54 


-124.94 


F8 


0.36 


0.17 


28.82 


40.30 


43.34 


2545.20 


F9 


0.36 


-0.17 


28.83 


40.30 


43.34 


-2545.20 


F10 


0.86 


0.00 


71.09 


38.95 


42.76 


0.00 


Fll 


0.53 


-0.33 


41.81 


40.07 


43.77 


-389.84 


F12 


0.53 


0.33 


41.81 


40.07 


43.77 


389.84 


F13 


0.20 


0.00 


18.59 


39.28 


43.02 


0.00 



Table 3. Changes in RMS WFE, image RMS width, and chro- 
maticity for the 6 selected cases discussed in sect. [5] 



Case 


RMS WFE 
[nm] 


Image RMS 
width [mas] 


Chromaticity 
[mas] 


1 


0.00 


39.26 


0.000 


2 


40.17 


42.51 


1.025 


3 


33.94 


41.75 


1.011 


4 


21.50 


40.14 


0.000 


5 


40.17 


42.23 


0.000 


6 


102.67 


360.72 


-0.799 



are based on improving the general image quality parameters 
such as WFE and spot diagrams. Therefore, it is possible to in- 
troduce custom merit functions in the optimisation procedure, 
which can include computation of the chromaticity or of the 
critical anti-symmetric contribution to WFE, with some aver- 
aging rule over the field. This does not necessarily modify the 
actual optical configuration in any significant way, since it is al- 
ways convenient to start after standard optimization, but chro- 
matic aberrations can be further reduced at the expense of the 
others and, possibly, of a small increase in overall WFE, which 
is acceptable in many cases. 

6. Field distribution of chromaticity 

The map of chromaticity over the FOV was derived from the 
optical representation of the baseline configuration of the astro- 
metric payload. We considered only the nominal design aberra- 
tions; manufacturing and alignment errors are not included yet. 
The result can therefore be considered as a best case. 
The analysis used monochromatic PSFs in the set of positions 
listed in Table |2 that also lists the associated RMS WFE in 
nanometers (nm), the image RMS width for the two cases con- 
sidered (B3 V and M8 V sources), and the chromaticity. The val- 
ues were interpolated to cover the field with 0°.02 resolution, 
and the resutling surface is shown in Fig. [6] The distribution is 
antisymmetric with respect to the y (across scan) axis, due to 
the intrinsic symmetry of the optical configuration: the mean 



T 




-40 -20 20 

along scan position (H m ) 



Fig. 5. Modifications of the PSF profile due to the different 
weight of aberration terms for field position F4 of the Gaia-like 
astrometric optics considered here. 



chromaticity value is -0.05 ^as, with RMS value 1.19 mas and 
peak values exceeding ±3 mas, i.e. three orders of magnitude 
larger than the Gaia measurement goal. 

6.1. Transit level compensation 

The symmetry of the chromaticity distribution can be exploited 
to reduce the overall contribution to a set of measurements, at 
transit-level, in spite of comparably high local values. All tar- 
gets detected by Gaia are observed in TDI mode on the whole 
FR so that the measurements are performed in opposite po- 
sitions with respect to the symmetry axis; by composition of 
the photo-center values from exposures in symmetric positions 
along scan, the residual chromaticity cumulated over a transit 
drops to values quite close to zero. 

The average value of the transit chromaticity over the astro- 
metric field is -0.10 fias, with RMS value 0.21 fias. Therefore, 
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Table 4. Transit-averaged chromaticity as a function of re- 
alignment error. 




Across [degrees] 



0.2 -0.4 



Along [degrees] 



Fig. 6. Chromaticity distribution over the astrometric field in 
the nominal configuration of the Gaia-like astrometric payload. 

chromaticity compensation over the transit appears quite effec- 
tive in the nominal configuration. 

6.2. Partial compensation on misaligned systems 

The real configuration used in flight will not retain the design 
symmetry, due to manufacturing, integration, and re-alignment 
errors. In order to model the effect, the chromaticity map de- 
rived for the nominal configuration is simply shifted by an 
amount corresponding to the selected offset ranging from one 
arcsecond to 0.5 arcminute (a fairly large alignment error). The 
local chromaticity value is mostly preserved, but the transit 
level combination, which is very close to zero in the nominal 
(symmetric) case, is affected by a degradation increasing with 
the error, as shown in Fig. [7] The statistics of transit-level chro- 
maticity across the FP is shown in Table |4] Column 1 lists the 
offset applied (zero corresponds to the nominal case); in cols. 2 
and 3 we report average and standard deviation of transit-level 
chromaticity, respectively, as computed across the field. 
The transit-level chromaticity remains very close to zero in 
an across-scan position of about 0°.55, which appears to be 
a chromatic-free section of the field; the residual has opposite 
signs on either side of this position. This is due to the struc- 
ture of the local chromaticity distribution (Fig. |6j, with alter- 
nate signs in each quadrant. The result is that stars of a given 
spectral type in different regions of the field have nearly op- 
posite residual chromaticity. This aspect may be exploited for 
further reduction of the residual chromaticity in the data re- 
duction phase. Upon definition of a threshold of acceptable 
residual chromaticity, at transit level, it is possible to provide 
a specification for alignment, referred to both initial telescope 
integration, and in-orbit re-alignment. 

6.3. Optical engineering aspects 

Symmetry considerations can be applied to the contribution 
from the individual optical components, in particular concern- 
ing their figuring error in the manufacturing phase, which can 
be split in components relevant or not to overall chromaticity. 
Also, alignment aspects, above referred to the primary mirror, 



Offset 
(arcsec) 


Chromaticity Mean 
(fias) 


Chromaticity RMS 
(juas) 





-0.403 


0.8218 


1 


-2.599 


14.1906 


5 


-11.277 


68.3465 


10 


-21.892 


136.5464 


30 


-61.998 


414.4432 
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Fig. 7. Transit level chromaticity vs. across scan field position 
in the nominal case and for increasing alignment error, (in arc- 
sec). 



can be applied with appropriate modifications to perturbation 
of other components. Manufacturing and alignment errors in- 
ducing an anti-symmetric residual WFE on a telescope mir- 
ror result in increased chromaticity. The on-flight configura- 
tion will be affected by manufacturing and alignment errors, 
which can also be described in terms of symmetric and anti- 
symmetric aberrations. If the design is optimized in terms of 
low chromaticity, small perturbations due to manufacturing and 
alignment are likely to induce a comparably low increase in 
the chromatic errors, due to larger tolerances. With the con- 
figuration considered as representative for Gaia, the relevant 
values of residual transit chromaticity are on the order of a few 
arcseconds at the level of the primary mirror, which is chal- 
lenging for standard optical techniques. However, it is not nec- 
essary to perform re-alignment of each optical component: as 
for many astronomical instruments, when the errors are not too 
large, a global re-optimization can be achieved by adjusting a 
small number of degrees of freedom, often localized on a single 
component (e.g. the secondary mirror). The correction require- 
ments in terms of stroke and resolution are within the range of 
the actuators considered for Gaia; the main limiting factor is 
the diagnostics capability, which in any case will benefit from 
accurate analysis of the image properties over the field. 
Quantitative analysis and optimization may be performed on 
specific optical configurations, including manufacturing and 
alignment aspects based on realistic WFE data from manu- 
facturers. The resulting minimization of the instrumental chro- 
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maticity is an improvement to the overall systematic error bud- 
get, which is desirable because correction procedures based on 
the science data are necessarily limited, and in some ways they 
subtract information. 

7. Conclusions 

This paper has discussed the issue of chromaticity in all- 
refiective optical systems. In spite of common belief, avoiding 
refractive components is not sufficient for achieving an achro- 
matic instrument response, due to basic diffraction consider- 
ations. In the scenario of future space astrometric missions, 
the impact of systematic errors two or three orders of magni- 
tude larger than the measurement goal is of fundamental im- 
portance. 

Specific analyses have been performed with reference to the 
baseline configuration of the Gaia astrometric payload, but the 
assumptions, principles, and conclusions of our discussion can 
be applied to any high accuracy astrometric instrument. The 
chromatic error is defined as the difference in image photo- 
center location at different wavelengths, and the exact value 
depends on the selected measurement process, but the effect is 
unavoidable. 

Independently from the selected WFE expansion, the terms rel- 
evant to chromaticity are those associated with anti-symmetry 
of the PSF in the FP and with anti-symmetric WFE contribu- 
tions on the pupil plane. Symmetric terms only contribute to 
the astrometric noise by increasing the effective image width in 
the measurement direction. As the relation between RMS WFE 
and chromaticity is complex, the specification of only the RMS 
WFE is not a sufficient requirement for controlling chromatic- 
ity from the optical manufacturing standpoint. 
The first prescription to optical manufacturers is to suppress 
or at least minimize the anti-symmetric terms. However, it is 
not necessary to set all chromatic terms to zero: an appropri- 
ate combination is still able to provide some local balancing. 
Standard optical design optimization techniques are able to 
provide reasonable results by applying the conventional image 
quality merit functions; optimal results on chromaticity require 
definition of ad hoc criteria, as much as on tolerancing. The dis- 
tribution of chromaticity over the field inherits some symmetry 
properties from the optical system; deviations from symmetry 
are induced e.g. by manufacturing and alignment errors on each 
optical component. 

In case of repeated measurements in different parts of the 
field, some chromaticity compensation is achieved in the data 
combination, depending on the symmetry of both instrument 
and measurement schemes. In the case of Gaia, each object 
observed by a symmetric (i.e. correctly re-aligned) telescope 
along a full transit provides a set of astrometric measurements 
affected by opposite chromatic errors in symmetric positions; 
transit-level composition is therefore likely to remove a large 
fraction of the local chromaticity. Any residual chromaticity 
must be removed in the science data processing after the best 
implementation of the astrometric payload, to minimize the ini- 
tial systematic error. This requires spectral information for each 
source and a good knowledge of the detailed instrument re- 
sponse (IGai & Cancellierel 2005). 
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